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ABSTRACT 

Possible  sources  of  the  marked  dicrepancy  between  empirical  and 
theoretical  (paticularly  electron-heated)  models  of  the  solar  chromosphere 
during  flares  are  discussed.  A  major  source  of  uncertainty  in  empirical 
models  is  the  inhomogeneity  of  the  spectral  data  on  which  they  are  based. 
Theoretical  models  involve  several  possible  sources  of  error  of  which 
probably  the  most  important  is  neglect  of  the  radiative  coupling  of  upper 
and  lower  chromospheric  regions.  In  addition,  both  types  of  model  are 
limited  by  their  assumption  of  plane-parallel,  and  frequently  static, 
structures. 

The  variety  of  deconvolutions  involved  in  the  construction  and  inter¬ 
comparison  of  empirical  and  theoretical  models  is  another  source  of  diffi¬ 
culty  in  model  testing.  A  new  procedure  for  studying  flare  energy  input 
is  suggested  whereby  the  required  input  is  derived  from  the  empirical 
model  chromosphere;  this  procedure  is  applied  to  the  electron-heated 
case,  and  it  is  found  that  the  integral  equation  defining  the  flare  energy 
deposition  rate  (erg  g  ^  s  S  can  be  analytically  inverted  to  yield  the 
injected  electron  flux  energy  spectrum  from  knowledge  of  the  energy  balance 
in  the  empirical  atmosphere.  We  analyze  recent  empirical  model  results  in 
this  manner  and  compare  the  calculated  injected  electron  flux  spectrum  with 
that  needed  for  hard  X-ray  bursts  in  moderately  large  flares.  This  technique 
betters  the  usual  procedure  of  predicting  a  model  chromosphere  from  the 
hard  X-ray  electron  flux  in  that  it  yields  a  measure  of  significance  and 
uniqueness  of  the  heating  model  parameters. 

Subject  Headings:  hvdromagnetics  -  Sun:  chromosphere  -  Sun:  flares  - 


Sun:  X-ravs 


I  INTRODUCTION 


A  considerable  effort  in  flare  physics  has  recently  gone  into  the  study 
of  the  secondary  response  of  the  solar  atmosphere  to  the  primary  flare  energy 
release  (see^.g.  ^Canfield  et  al.  1980  for  a  review) .Essentially  this  field 
comprises  the  comparison  of  flare  model  atmospheres  constructed  in  two  dis¬ 
tinct  ways,  viz.  theoretically  and  empirically.  A  theoretical  model  is  an 
atmospheric  structure  predicted  by  solution  of  the  equations  (e.g.  transport 
of  energy,  momentum  and  radiation)  describing  the  response  of  the  preflare 
solar  plasma  to  a  hypothetical  flare  energy  input.  An  empirical  model,  on 
the  other  hand,  has  a  structure  derived  from  observations  of  the  flaring 
atmosphere,  by  trial  and  error  fitting  of  the  computed  model  spectrum  and 
subsequent  model  adjustment.  The  theoretical  model  may  then  be  compared  with 
the  empirical  or  (which  is  not  entirely  equivalent  -  see  §111)  be  used  to 
predict  a  spectrum  for  direct  comparison  with  observations.  Recent  work 
(Machado  and  Emslie  1979  ;  Lites  and  Cook  1979  ;  Machado  et  al.  1980 — hereafter 
MAVN)  on  empirical  models  has  established  marked  discrepancies  between  these  and 
theoretical  models,  particularly  those  inwhich  the  hypothetical  energy  input 
is  particle  beam  collisions  (e.g.  Brown  1973,  Brown,  Canfield  and  Robertson 
1978 — hereafter  BCR;  LaBonte  1978;  Machado,  Emslie,  and  Brown  1978).  Not  only 
is  there  a  quantitative  discrepancy  in  the  parameters  required  to  describe  the 
theoretical  energy  input  (e.g.  electron  flux)  and  the  observations  (e.g. 

Ha  brightness)  of  part  of  a  particular  flare  (see  BCR)  but  also  qualita¬ 
tive  differences  between  the  general  form  of  empirical  and  theoretical 
structures  for  any  input  parameters.  In  this  paper  we  discuss  possible 
sources  of  these  discrepancies  with  a  view  to  clarifying  the  most  important 
areas  for  future  work  on  this  problem.  Our  discussion  is  concerned  mainly 
with  models  of  the  chromosphere,  and  specifically  electron-heated  models. 
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The  problem  of  the  optically  thin  corona  and  transition  region  has  been 
discussed  recently  by  Machado  and  Emslie  (1979)  and  that  of  the  tempera¬ 
ture  minimum  region  by  Machado,  Emslie,  and  Brown  (1978).  However,  many 
of  our  remarks  on  the  sources  of  discrepancy  in  the  chromosphere  are 
equally  applicable  to  these  regions  and  to  models  heated  by  other  forms 
of  energy  input,  such  as  soft  X-ray  irradiation  (Somov  1975;  Machado  1978). 

Figure  1 

In  Figure  1  (from  MAVN  )  we  compare  the  empirical  models 

derived  from  Sky lab  ATM  data  with  the  electron  heated  (BCR)  theoretical  models. 
The  main  discrepancies  are: 

(i)  The  BCR  models  have  higher  temperatures  in  the  mid-chromosphere 

(T  ^  104K,  h  -  1100  km)  and  much  more  material  with  104  <  T  <  2  x  104K 

r\j 

(ii)  The  transition  region  in  the  BCR  models  occurs  considerably  higher 
in  the  atmosphere  than  in  the  empirical  models.  Indeed  for  small 
electron  fluxes  the  BCR  transition  region  occurs  higher  in  the 
flaring  atmosphere  than  in  the  quiet  sun  (Machado  and  Emslie  1979). 

In  the  region  h  -  800-1100  km  (T  3  8  *  103  -  104K)  the  BCR  and  empirical 
models  have  similar  temperatures;  to  that  extent  the  BCR  models  are  locally 
plausible,  reproducing  observed  Ha  profiles  when  a  moderate  amount  of 
horizontal  inhomogeneity  (see  §IIc),  consistent  with  the  filling  factor  of  the 
Ha  spectrograph  slit,  is  included.  (Note  however  that  the  electron  flux 
needed  to  obtain  this  agreement  is  ^  1/10  of  that  required  for  thick  target 
interpretation  of  hard  X-rays  from  the  same  flare  -  BCR.)  Nevertheless  such 
a  spatially  localized  regime  as  that  of  Ha  formation  clearly  cannot  be 
expected  to  provide  an  adequate  empirical  test  of  a  model  energy  input. 
Furthermore  (see  §IIb)  a  localized  treatment  is  not  satisfactory  theor¬ 
etically  either.  It  is  therefore  essential  to  take  a  global  view  in 
modeling  of  the  chromosphere,  and  indeed  of  the  whole  atmosphere. 

In  5 1 1  we  consider  the  possible  sources  of  the  discrepancies  outlined 
above  and  discuss  the  global  radiative  transfer  problem  more  fullv.  In 
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§  III  we  discuss  the  nature  of  the  intercomparison  process  itself 
and  how  it  renders  estimation  of  uncertainties  difficult.  We  then 
suggest  a  new  intercomparison  technique  which  allows  estimation  of  these 
uncertainties;  this  technique  involves  computing  the  required  energy 
input  for  a  given  empirical  atmosphere  structure,  rather  than  constructing 
a  theoretical  model  atmosphere  based  on  a  prescribed  hypothetical  flare 
energy  input.  In  §IV  we  apply  this  technique  to  MAVN  empirical  flare 
data  and  discuss  the  results  obtained. 


-fe¬ 
ll.  SOURCES  OF  DISCREPANCY 

a)  General 

There  are  several  features  in  common  to  both  empirical  and  theoretical 
models  which  render  them  generally  unsatisfactory.  First ,  in  order  to  make 
the  analyses  tractable,  most  chromospheric  modeling  has  thus  far  been  based 
on  the  assumption  of  a  steady  state  plane-parallel  structure  in  hydrostatic 
equilibrium.  These  assumptions  eliminate  time-dependence  from  the  equations 
and  permit  representation  of  the  models  purely  in  terms  of  their  temperature 
versus  height  structure  T(h)  as  used  in  Figure  1,  together  with  a 
scaling  factor  defining  pressure  or  density  at  a  single  point.  However, 
they  result  in  omission  of  energy  transport  and  spectroscopic  effects 
due  to  mass  motion,  time  dependent  conduction  (see  Craig  and  McClymont 
1976)  and  radiative  transport  (see  Dubov  1963),  and  inhomogeneity.  For 
example,  if  large  velocity  gradients  were  present  in  the  atmosphere, 
then  it  is  possible  that  they  would  radically  alter  the  interpretation 
of  spectral  diagnostics  in  terms  of  a  T(h)  structure.  (There  is,  however, 
little  evidence  for  such  velocity  gradients  in  the  layers  of  the  flaring 
atmosphere  under  consideration  [see  Lites  and  Cook  1979}.)  Second, 
chromospheric  modeling  of  both  types  involves  the  complex  numerical  solution 
of  the  relevant  radiative  transfer  equations,  the  techniques  of  which  are 
still  under  development  (see  Vernazza,  Avrett,  and  Loeser  [VAL]  1973, 

1976,  1980)  and  thus  require  cautious  treatment. 

b)  Deficiencies  in  Theoretical  Models 


A  number  of  factors  in  the  electron  beam  energy  input  used  by  BCR  warrant 


further  study.  First,  the  input  is  purely  collisional  (see  Brown  1973; 
Emslie  1978)  taking  no  account  of  the  effects  of  plasma  wave  generation 
i Brown  and  Melrose  1977;  Vlahos  and  Papadopoulos  1979)  nor  of  reverse 
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current  ohmic  losses  (e.g.  Emslie  1980)  on  the  height  distribution  of 
beam  energy  deposition.  These  are  both  most  important  highest  in  the 
atmosphere  and  will  tend  to  steepen  the  BCR  T(h)  profiles  toward  the 
empirical  forms,  both  directly  and  by  enhancing  the  role  of  conduction 
(see  below).  Second,  the  effect  on  beam  dynamics  of  converging  magnetic 
field  lines  is  neglected;  its  inclusion  would  enhance  the  beam  deposition 
at  higher  altitudes.  Third,  the  BCR  models  were  computed  only  for  electron 
spectral  index  <S  *  4  .  Use  of  the  higher  5's  typical  of  most  small 
flares  will  also  tend  to  steepen  T(h)  by  depositing  energy  preferentially 
higher  in  the  atmosphere  (see  Emslie  1978). 

The  BCR  models  included  no  thermal  conduction  on  the  grounds  that 
conductive  heating  (or  cooling)  is  negligible  at  the  low  temperatures 
and  temperature  gradients  where  Ha  forms  (see  Brown  1974,  1977).  We 
have  checked  this  assumption  directly  for  both  the  theoretical  and  empiri¬ 
cal  models  and  find  (see  Sill)  that  for  T  <  12000  K  the  conductive  term 


dh 


<(T) 


dT 

dh 


is  totally  negligible  compared  to  the  radiative  losses  there  so  that  an 
external  input  is  needed.  This  is  contrary  to  claims  by  Svestka  (1970), 
Shmeleva  and  Syrovatskii  (1973),  and  Davis,  Kepple,  and  Strickland  (1977) 
that  local  conductive  deposition  can  balance  chromospheric  radiative 
losses  during  the  post-flash  phase  of  the  flare.  At  first  sight  this 
also  appears  contrary  to  the  statement  by  Machado  and  Emslie  (1979)  that 
"thermal  conduction  alone  is  sufficient  to  produce  the  observed  flare 
enhancements  at  all  atmospheric  levels,  in  particular  in  the  chromosphere". 
However,  it  turns  out  that  in  the  proper,  global,  approach  to  the  theoreti¬ 
cal  modeling,  it  is  essential  to  include  conduction  in  the  hotter  layers 


(upper  chromosphere)  in  order  for  the  correct  structure  in  the  lower  (Ha) 
chromosphere  to  be  obtained.  The  reason  is  that  these  two  regions  are 
radiatively  coupled,  particularly  by  La  which  constitutes  an  important 
energy  loss  from  the  upper  chromosphere  but  a  significant  input  into  the 
lower  chromosphere  —  Machado  and  Emslie  (1979);  MAVN.  (The  possibility 
of  radiative  heating  was  noted  by  BCR  but  not  incorporated  in  their 
models.)  The  importance  of  this  La  "backwarming"  stems  from  the  fact 
that  the  energy  in  the  La  radiation  comes  from  dissipation  of  thermal 
conductive  flux,  where  the  steep  transition  zone  temperature  gradient  is 
suddenly  reduced.  In  short,  the  Machado  and  Emslie  (1979)  statement  refers 
to  the  fact  that  a  thermal  conductive  flux  in  the  transition  zone  can, 
through  La  backwarming,  supply  the  energy  for  a  significant  fraction 
of  the  radiative  losses  in  the  upper  chromosphere  (i.e.  down  to  about  the 
region  where  Ha  forms) . 

The  excessively  high  BCR  transition  region  is  attributable  substantially 

to  two  factors.  First  BCR  take  the  quiet  sun  model  of  VAL  1973  as  their 

preflare  state  in  establishing  the  density  scale  and  the  ambient  energy  input 

in  their  models  (see  Brown  1973).  Use  of  a  plage  preflare  model  would  have 

resulted  in  deeper  transition  regions,  particularly  in  low  flux  models.  The 

residual  discrepancy,  and  in  particular  the  fact  that  the  BCR  transition  zone 
for  low  electron  fluxes  lies  above  its  quiet  sun  level,  is  a  numerical 

artifact.  Since  BCR  were  attempting  to  model  the  Ha  formation  region,  and 

since  they  did  not  recognize  the  importance  of  the  radiative  coupling  of 

upper  and  lower  chromospheric  regions,  they  used  only  a  very  broad  grid 

representation  of  the  upper  chromosphere.  In  fact  the  grid  step  used 

there  was  of  the  order  of  the  discr  -.ncy  in  transition  zone  height  noted 


by  Machado  and  Emslie  (1979). 


c)  Deficiencies  in  Empirical  Models 

Aside  from  the  points  raised  in  §IIa,  and  the  problems  of  model  unique¬ 
ness  discussed  below  (illl),  the  main  deficiencies  in  empirical  models  lie  in 
the  inhomogeneity  of  the  data  used.  This  takes  several  forms.  First,  the 
evident  inhomogeneity  of  the  horizontal  structure  of  the  flare  within  the 
instrumental  field  means  that  the  empirical  models  represent  some  horizontally 
averaged  flare  structure.  It  is  therefore  not  particularly  meaningful  to 
compare  these  with  theoretical  models  representing  the  flare-heated  areas 
(as  MAVN  do  in  Figure  1).  Rather  we  should  compare  them  with 

a  theoretical  model  diluted  by  the  ambient  atmosphere  with  an  appropriate 
filling  factor  (which  itself  generally  depends  on  height  h  according  to 
the  variation  of  flare  area  and  of  instrumental  field  with  h  )  .  The 

latter  is  the  procedure  followed  by  BCR  in  their  Ha  analysis. 

Second,  even  the  best  observational  efforts  to  date  (  MAVN  ) 
have  resulted  in  spectra  based  on  very  inconsistent  coverage  in  both 

space  and  time.  In  some  cases  the  data  used  in  different  spectral  bands 
were  not  even  from  a  single  flare,  and  at  best  were  from  different  parts  of 
a  single  flare  at  different  times.  Furthermore  much  of  the  ATM  flare  cover¬ 
age  was  in  the  post-flash  phase,  so  putting  in  doubt  the  relevance  of  the 
empirical  models  in  testing  theoretical  models  heated  by  electron  beams 
which  are  early-stage  phenomena. 

Third,  the  flare  set  analyzed  by  MAVN  comprised  five 

subflares  and  one  Class  L,  while  some  of  the  theoretical  model  analyses 
(e.g.  BCR)  have  been  aimed  at  very  large  events. 

III.  CONVOLUTIONS  AND  STABILITY  IN  MODEL  BUILDING 
An  unsatisfactory  aspect  of  current  empirical  models  is  that  they  give  no 


error  bars  on  T(h)  so  that  any  comment  on  their  compatibility  with  theore- 
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tical  models  is  only  qualitative.  The  question  we  really  want  to  ask  is:  what 
manifold  of  empirical  model  atmospheres  is  capable  of  fitting  the  spectral 
data  within  their  error  bars,  and  what  theoretical  models  fall  within  this 
manifold?  The  first  part  of  this  question  is  hard  to  answer  due  to  the  very 
large  computational  effort  involved  in  iteration  to  obtain  even  one  acceptable 
empirical  model.  Furthermore  error  bars  on  T(h)  so  obtained  are  almost 
certain  to  be  rendered  large  and  hard  to  interpret  meaningfully  because  of 
the  inhomogeneity  of  data  used.  Also,  in  comparing  theoretical  and  empirical 
models  quantitatively ,  we  should  develop  an  objective  criterion  of  acceptabi¬ 
lity  .presumably  based  on  some  statistic  such  as  the  root  mean  square  residual 
on  the  energy  balance  parameters  (rather  than  simply  on  temperature)  relative 
to  the  probable  error  in  the  observed  (empirical)  values. 

Even  given  a  good  set  of  homogeneous  data  with  error  bars,  there  are 
still  fundamental  limitations  on  the  extent  to  which  empirical  models  can 
test  theoretical  models  because  of  the  number  of  deconvolution  processes 
involved  (see  Craig  and  Brown  1976)  .  Specifically,  in  testing  electron- 
heated  chromospheric  models,  the  usual  procedure  is  (a)  derive  an  electron 
spectrum  from  a  noisy  observed  hard  X-ray  spectrum  by  deconvolution  (Brown 
1971;  Craig  1979);  (b)  predict  a  theoretical  model  atmosphere  heated  by 

this  electron  input  by  solving  the  integro-dif ferential  equation  of  radiative 
transport;  (c)  derive  an  empirical  model  atmosphere  structure  by  decon¬ 
volution  of  the  integro-dif ferential  equation  governing  the  formation  of 
the  observed  spectrum;  (d)  compare  (c)  with  (b)  .  Given  the  problems 
of  error  magnification  in  deconvolut ing  from  and  of  nonuniqueness  in 
model-fitting  to,  integral  equations  (Craig  and  Brown  1976  and  references 
therein),  we  should  clearly  be  circumspect  in  rejecting  or  accepting 
theoretical  model  chromospheres  heated  by  electrons,  since  their  computa¬ 
tion  depends  intrinsically  on  the  solution  of  t hree  integral  equations. 


Here  we  suggest  the  following  alternative  to  the  usual  sequence  above, 
which  enables  us  to  quantify  uncertainties  in  derived  quantities; 
the  usual  process  of  forward  fitting  of  models  to  data  does  not  permit 
us  to  do  this. 

Instead  of  predicting  an  optical/UV  spectrum  from  a  theoretical 
model  flare  chromosphere  with  the  energy  input  function  inferred  from 
independent  observations  and  comparing  this  prediction  with  optical  and 
UV  observations,  we  suggest  reversing  the  procedure.  That  is,  we  derive 
the  energy  input  function  which  is  consistent  with  the  heating  rates 
derived  from  the  optical/UV  spectral  data  and  compare  this  with  other 
independent  estimates  of  this  input  function.  (See  Machado,  Emslie, 
and  Brown  1978  for  similar  discussion  of  heating  at  the  temperature  minimum.) 
Specifically  we  propose  the  following  procedure: 

(a)  Obtain  empirically  that  atmospheric  temperature  structure  T(h) 
which  best  fits  the  spectral  data,  or  better  that  range  of  T(h) 
which  fit  the  data  within  its  errors. 

(b)  Evaluate  the  net  radiative  energy  loss  from  this  (these)  empirical 

model (s)  (see  VAL  1980)  and  so  deduce  the  total  specific  energy 

-1  -1 

input  I  (h)  (erg  g  s  )  required  for  a  steadv  state.  (I  [h ] 

°  '  o 

comprises  a  gross  radiative  loss  I  th ]  less  a  radiative  input 

tot  r 

I rf h ]  from  backwarming. ) 

(c)  Subtract  from  I  (h)  the  specific  conductive  input  I  *t  !  X  i-  (vrrl  — ) 

°  c  (pdhllldh. 

(p  =  gas  density)  calculated  from  the  T(h)  found  in  (a),  and  also  an 
estimate  of  the  ambient  energy  input  Ia(h).  The  latter  requires  some 
assumption  of  how  Ia(h)  changes  between  preflare  and  flare  condi¬ 
tions — e.g.  constant  I  per  gram  (see  Brown  1973). 

d 
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(d)  This  leaves  the  energy  input  I^(h)  to  be  supplied  by  the  flare 
release  processes,  which  can  then  be  compared  with  theoretical 
heating  models.  The  acceptability  of  the  fit  as  the  various  para¬ 
meters  of  the  theoretical  model  are  varied  gives  an  indication  of 
the  uncertainties  in  the  values  of  these  parameters,  quantities  which 
are  difficult  to  obtain  in  the  conventional  comparison  procedure 
( see  BCR) .  It  should  also  be  noted  that  this  procedure  should 
enable  a  direct  comparison  of  1^  with  1^,  I^_  and  1^  at  each 
height  h  ,  indicating  unambiguously  the  relative  role  of  each 
process  (as  a  function  of  height),  which  up  to  now  has  remained 
controversial . 

We  illustrate  the  technique  by  considering  chromospheric  models  heated  by 
collisional  degradation  of  a  beam  of  non-thermal  electrons,  injected  vertically 

along  uniform  field  lines.  Expressing  1^  as  a  function  of  the  Lagrangian  column 
depth  variable  N  =  r  n(h')dh'  rather  than  h  ,  then  the  energy  input 


arising  from  an 


injected  number  flux  spectrum  F  (E  )  (electrons 

o  o 


cm  2 s  1  per  unit  Eq)  is  (Brown  1973;  Ernslie  1978) 


F  (E  ) dE 


^1  J  /3KN  Eo(1"3 


KN/E2) 2/3 


(1) 


where  K  =  2Tte  A  in  the  usual  notation  and  m^  is  the  mass  of  a  hydrogen 

K  F  (/x) 

atom.  Setting  x  =  ,  v  =  3KN,  G(x)  =  - and  J(v)  *  I.(N)  then 

H  x*  1 


J(,) . /-  , 


(2) 
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which  is  an  Abel  integral  equation  with  solution  (assuming  G(x  - 
all  reasonable  F  (E  )) 


G(x) 


[4 


J(v)  +  vJ*  ( 


-1 


dv 


(v-x) 


1/3 


0  for 


(3) 


(  Courant  and  Hilbert  1953,  pp.  158—159) - 
In  terms  of  the  original  variables* 


F 

o 


(Eo)  =  - 


2  ^  r- 

wE2  *42/3K 

o  o 


|(If(N)  +  \  N  If(N)JdN 


(f-r 


(4) 


which  determines  the  form  of  the  injected  electron  flux 

spectrum  Fo(Eq)  from  knowledge  of  the  behavior  of  1^  with  N,  obtained  via 
steps  (a)  to  (c)  above.  In  the  next  Section  we  shall  apply  this  method  of 
analysis  to  empirical  flare  data  from  MAVN. 


IV  APPLICATION  TO  OBSERVATIONS 


In  Figure  2  we  show  the  behavior  of  I  ,  I  ,  I  and  I  obtained  by 

o  tot  r  c 

the  procedure  of  VAL  (1980)  for  model  MAVN2  (see  Figure  1) .  We  use  this 
model  (as  opposed  to  model  MAVN1)  because  it  corresponds  to  bright  flare 
regions  and  so  is  more  comparable  with  an  atmosphere  undergoing  excitation 
by  non-thermal  particles,  though  it  is  still  small  compared  to  a  BCR-sized 
flare.  One  immediately  notes  the  following  points  from  the  Figure: 

(i)  the  local  condutive  input  Ic  is  less  than  Iq  by  at  least  nine 
orders  of  magnitude  and  so  is  totally  negligible  (see  remarks  in 
§11). 

(ii)  the  ratio  I  /I  (the  contribution  of  backwarming  to  the 

r  tot 

radiative  energy  budget)  approaches  unity  (to  within  a  factor 

of  -  2;  see  remarks  in  §IIb)  at  around  N  -  3X1021  cm  2,  but 

falls  off  rapidly  with  depth  thereafter  (so  that  I  -  I  for 

tot  o 

depths  greater  than  this). 

To  derive  I^(N)  we  further  need  to  apply  process  (c),  i.e.  subtract 
the  ambient  (preflare)  input  I  (N)  and  the  local  conductive  input  I  (N) 
As  noted  in  (i)  above,  the  latter  of  these  is  negligible.  1^  can  be 
obtained  by  applying  the  VAL  (1980)  method  of  calculating  net  radiative 
losses  to  a  plage  model  atmosphere,  preferably  an  empirical  one.  The  only 
available  empirical  plage  model  is  that  due  to  Basri  et  al.  (1979);  however 
this  model  suffers  from  serious  errors  in  the  upper  chromospheric  region, 
due  to  the  assumption  of  local  thermodynamic  equilibrium  for 

the  Carbon  atom  (E.  H.  Avrett,  private  communication).  Since  this  region 
is  precisely  the  region  of  greatest  interest  to  our  present  study,  we  do 
not  consider  it  appropriate  to  use  the  Basri  et  al .  (1979)  model  to  obtain 
I  (N’)  .  We  therefore  have  no  alternative  but  to  neglect  I  (N)  and  res- 

3  3  ~ 
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Figure  3 


l 


reasonably  taken  to  be  negligible  compared  to  I  ;  hopefully  future  work  on 
empirical  plage  model  atmospheres  may  enable  us  to  remove  this  restrictive 
condition.  Figure  2  strongly  suggests  that  I  becomes  important  around 
N  -  1022  cm  2  (note  the  upturn  in  the  graph  of  Iq(N)  at  this  point);  this 
belief  is  somewhat  strengthened  by  application  of  the  VAL  (1980)  radiative 
loss  calculation  to  the  Basri  et  al.  (1979)  model,  which  shows  that  radiative 
losses  due  to  Mg  and  Ca  lines  start  to  become  very  important  at  just  this 
depth. 

In  Figure  3  we  show  the  result  of  applying  equation  (4)  to  the  Iq 
(=  I  ,  since  we  are  neglecting  I  )  values  of  Figure  2.  to  avoid  the 

I  3 

singularity  at  the  lower  integral  limit  in  our  computations  we  first  trans¬ 


formed  (4)  to  the  equivalent  expression 


F  (E  ) 
o  o 


[l+x3/2] 


) 


dx 


(5) 


(see  Brown  1971) .  In  the  evaluation  of  this  integral  there  is  a  truncation 

error  at  the  upper  limit,  since  we  only  have  data  up  to  x  *  (3KN  /E2-l)2/3 
r  max  max  o 

where  N  is  the  largest  column  density  for  which  data  is  available.  This 
max  ® 

truncation  error  will  only  seriously  affect  the  determination  of  F  (E  )  near 

o  o 

E  =  ( 3KN  )h  ,  which  is  well  outside  the  range  of  E  values  for  which  a 
o  max  o 

spectrum  is  to  be  determined  (see  below);  thus  it  does  not  affect  the  results 

to  follow.  Note  the  hardening  in  the  electron  spectrum  above  around  300  keV, 

22  -2 

corresponding  to  N = 10  cm  —  cf.  the  lower  integral  limit  in  equation  (4). 

This  is  attributable  in  part  to  our  neglect  of  I  ,  as  noted  above,  but  is 

a 

certainly  also  in  part  representative  of  the  strong  heating  observed  near  the 
temperature  minimum  which  cannot  be  explained  by  non-thermal  electron  bombard¬ 
ment  (Machado,  Emslie,  and  Brown  1978).  Shown  in  the  Figure  is  the  best-fit 
power-law  F  (E  )  *  F  E  ~6  to  the  (E  ,F  )  points,  where  only  E  values  to  the 

O  O  OO  O  O  O  '  o 

left  of  the  spectral  bend  have  been  used  (see  remarks  above  on  the  limits  of 
applicability  of  our  method).  This  best-fit  line  has  the  form 


..A 
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F  (E  ) 
o  o 


10 


(26.91i  0.53)  --(6.21  ♦  0.24) 


—  2  —1 

elec crons  cm  s  keV 


(6) 


with  a  regression  coefficient  r  -  0.997  . 

For  comparison  with  hard  X-ray  burst  analysis  it  is  more  convenient  to 
express  F^o  in  terms  of  the  total  injected  energy  flux  above  20  keV  (say), 
viz. 

F  (erg  cm'2s'’ )  -  1.6  *  10“’/  F  (E  )  E  dE  .  (7) 

20  «/  o  o  o  o 

2  0 


Using  equation  (6)  in  equation  (7)  we  obtain,  for  the  injected  electron  flux 
parameters : 


log  F  -  12.01  t  o.53 
5  -  6.21  t  0.24 


(8) 


From  these  results  and  their  uncertainties  we  conclude  that  heating 
by  an  electron  beam  with  a  power  law  injection  spectrum  fits  the  observa¬ 
tions  extremely  well  over  the  atmospheric  domain  concerned;  however,  the 
scatter  in  the  few  data  points  used  prevents  us  from  deriving  the  required 
electron  flux  F,„  with  an  accuracy  better  than  about  half  an  order  of 
magnitude.  We  also  note  chat  the  value  of  5  in  eq .  (3)  is  much  larger 
than  the  •,»4  used  by  BCR.  This  explains  the  higher  temperatures  (compared 
to  the  MAVN  models)  found  by  BCR  in  the  Ha  formation  region  (il);  see 
comments  in  5 1 1  on  the  effect  of  different  5  on  the  variation  of  the 
energy  deposition  rate  with  depth.  With  this  in  mind,  we  note,  however, 
that  the  electron  injection  parameters  (8)  are  entirely  consistent  with 

those  required  for  hard  X-ray  bursts  in  large  events  (Hoyng,  Brown,  and 

-2-1 

van  Beek  1976),  viz.  F.  -  3  '  10‘*  erg  cm  s  (for  plausible  injection 
areas  A  -  10‘-  -  10‘-  cm-)  and  5-4-6.  This  is  encouraging,  although 
it  must  be  noted  that  the  above  illustrative  analysis  is  based  on  a  very 


limited  range  of  points,  viz.  a  factor  of  '-5  in  the  column  density  N 
and  so  a  factor  of  ^2  in  the  electron  energy  Eq  .  In  addition,  the 
error  on  F,,  quoted  in  eq .  (8)  is  based  on  the  residuals  of  the 
<Eo.Fq)  points  about  the  best  fit  line  (Figure  3)  and  might  be  substantially 
larger  if  the  uncertainties  in  the  individual  points  exceed  these  residuals; 
these  uncertainties  are,  however,  very  difficult  to  estimate  (E.H.  Avrett, 
private  communication) .  It  is  to  be  hoped  that  more  work  on  preflare 
structures,  and  the  acquisition  of  more  extensive  data  sets,  will  enable 
the  present  method  to  be  used  as  a  better  flare  diagnostic  technique. 
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FIGURE  CAPTIONS 

Figure  1:  Temperature  versus  height  structures  in  empirical  model  chromospheres 
for  small  (MAVN  1)  and  larger  (MAVN  2)  flares  constructed  from 
Skylab  data  by  Machado  et  al.  (1980),  compared  with  the  structures 
predicted  by  Brown,  Canfield,  and  Robertson  (1978)  for  chromospheres 

heated  by  electron  beams  with  power  law  spectral  index  4  and  energy 

9  10  -2-1 

fluxes  (above  20  keV)  of  10  (BCR  9)  and  10  (BCR  10)  ergs  cm  s 

respectively.  The  column  density  values  at  the  top  of  the  Figure 

correspond  to  model  MAVN  2  (cf.  Figure  2). 

Figure  2:  Energy  sources  and  sinks  as  functions  of  column  density  N 

inferred  from  the  empirical  MAVN 2  flare  structure.  I  (N)  is  the 

o 

net  radiative  loss  (ergs  g  s  ),  comprising  a  gross  radiative 
loss  term  I  and  a  (negative)  backwarming  term  I  .  I  (N) 

tOL  l  C 

is  the  conductive  deposition  inferred  from  the  temperature 
gradients  in  the  atmosphere.  Note  the  importance  of  backwarming 
effects  for  N  102z  cm  Z  and  the  negligible  amount  of  local 
conductive  heating. 

Figure  3:  Injected  electron  flux  spectrum  obtained  bv  inversion  of  the  I  (N) 

o 

values  in  Figure  2  (see  equation  [4]).  Beyond  Eq  =  300  keV 
(corresponding— see  equation  [4]  — to  N  =  10  cm~“)  the  spectrum 

changes  form  dramatically;  this  is  partly  real  and  partly  due 
to  the  neglect  of  preflare  heating  terms  in  the  flare  energy 
equation— see  discussion  in  text  (§111).  Below  this  break  point, 
the  spectrum  is  well  fit  by  a  power  law  with  the  parameters  shown; 
these  can  be  compared  to  independently  derived  electron  flux 
parameters,  obtained  from  analysis  of  hard  X-ray  spectra  (see  * 

equation  [8]). 

_  J 
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